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Abstract Stationary photon-atom entanglement is discussed by applying the 
method of flow equation to the Jaynes-Cummings model. A nonlocal continuous unitary 
transformation is explicitly constructed and the associated positive operator-valued 
measures for the photons and atom are obtained. Then, flow of the entanglement 
entropy is analyzed. A comment is also made on implementing the unitary operation in 
the method of flow equation. This method may offer a new strategy for quantum-state 
engineering. 
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I. INTRODUCTION 

A recent major breakthrough in experimental quantum information would be the 
realization of quantum teleportation between photons and atoms [1]. A crucial step in 
the experiment is to generate entanglement between polarization of light and caesium 
atoms by the use of the quantum-mechanical Faraday rotation. This development is of 
fundamental importance toward fault-tolerant quantum computation and quantum 
communication. 

Importance of photon-atom entanglement sheds fresh light on traditional quantum 
optical systems. Actually, there are some discussions about photon-atom entanglement 
in the literature. They are, however, concerned with nonstationary states with rapid 
temporal oscillations of entanglement [2,3], which may not be so useful for quantum 
information processing. Thus, a point is how to generate stationary entangled states in a 
photon-atom system. Mathematically, this may be nothing else than diagonalizing the 
system Hamiltonian. Experimentally, however, given an initial state, it is hard to realize 
the quantum operation that can immediately transform such a state of the coupled 
system to the exact stationary state, in general. Therefore, in reality, it is necessary to 
inquire into a series of "trial" operations in order to approach the stationary state. Then, 
the question is what a desirable strategy for it is. 

To answer to the above question, here we examine the method of flow equation 
introduced by Wegner [4,5]. This theory uses a continuous unitary transformation 
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depending on a single parameter, /. The Hamiltonian tends to become diagonalized in 
the limit, l^oo. It is our opinion that Wegner's flow equation has physical meanings 
more than just a procedure of diagonalizing the Hamiltonian. In fact, it may offer the 
optimal strategy for quantum-state engineering [6]. In other words, if the trial 
operations are ejficient, they are almost on a trajectory generated by the flow equation. 
Due to this observation, it seems meaningful to examine the flow equation approach for 
quantum optical systems of photons and atoms. 

In this paper, we therefore consider the flow equation approach to establishing 
stationary entanglement between light and matter. As a simple but nontrivial example, 
we study the celebrated Jaynes-Cummings model [7,8] and perform its complete flow- 
equation analysis. Then, we shall calculate the positive operator-valued measures 
(POVMs) for the photons and atom, evaluate the degree of stationarity, and analyze 
flow of the entanglement entropy. We also present a procedure for implementing the 
unitary operation appearing in the method of flow equation for the model. 

II. BASICS OF FLOW EQUATION 

Let us recapitulate the method of flow equation. Consider the continuous sequence of 
unitary transformations parametrized by /: U{1). I is referred to as the flow parameter, 
which has dimension (energy)"^. Application of U{1) to a Hamiltonian H yields 
H(l) = U (l)HU\l), which obeys the following differential equation: 
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al 



(1) 



with the "initial" condition H{Q) = H. ri(l) = [dU(l) I dl\ U\l) is an anti-Hermitian 
operator. Wegner's proposal for the choice of this operator is 

J](l) = \.H\l),H(l)], (2) 

where H'^(l) stands for the diagonal part of H(l). With this choice, it can be shown 
[4,5] that the off-diagonal matrix elements of H(l) exponentially decay in / if H'^(l) is 
nondegenerate. Eq. (1) defines the flow equations of the coefficients contained in H(l). 
This is in analogy with the concept of renormaUzation group. 

III. FLOW EQUATION OF THE JAYNES-CUMMINGS MODEL 

Let us apply the method of flow equation to the Jaynes-Cummings model [7,8]. This 
model describes quantum dynamics of the interaction between a single-mode radiation 
field and a two-level atom in the rotating wave approximation. The Hamiltonian reads 

H = —(j^ + (Da^a + A.((J^a + a_a^). (3) 
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The symbols appearing in this equation are defined as follows. = |e)(e|-|g)(g|, 
cr+ = \e){g\, o_ = \g){e\ are the operators of the two-level atom with the normalized 
ground and excited states, |g) and \e). and a are the creation and annihilation 
operators of the photon satisfying the algebra: [a,a^] = l, [a,a] = [a^,a^] = 0. COq, 
(0, and A are the frequency of transition between the ground and excited states of the 
atom, the frequency of the single-mode photon, and the real coupling constant between 
the photon and the atom, respectively, h is set equal to unity for the sake of simplicity. 
In the present work, we limit our analysis to the off-resonant case. Therefore, the 
detuning factor, A = ft),) - CO, is nonvanishing, and here it is assumed to be positive: 

A = «o-«>0. (4) 

Although this is a highly simplified model of the coupled system of light and matter, 
remarkably it can experimentally be realized to generate single-atom maser or 
micromaser in a cavity [9]. Since the model is exactly solvable, it allows to analytically 
discuss a variety of phenomena such as the Rabi oscillation and quantum 
collapses/revivals. 

For the Jaynes-Cummings Hamiltonian in Eq. (3), we have constructed the following 
operator of the continuous unitary transformation: 
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= [«„(0(7, G_\n){n\ + P„_,(l)(7_ (7,\n){n\ 

+ Y„(l)cj,\n){n + M + S„m_ \n + l){n\]. (5) 

Here, jS_j(/) = e'^*'^ [which turns out to be unity, later: see Eq. (27)] and 
|n) = (n\y^'^ (a^)"\0) denotes the n-photon state with the normalized vacuum state |0), 
i.e., a|0) = 0. The unitarity condition, U{l)U\l) = U\l)U{l) = I, leads to the 
following set of equations ( n = 0, 1, 2, • • 



i«„(or+ir„(or=iiS„(/)i'+i5„(/)r=i, 



|«„(0r + |<5„(/)|' = |/?„(/)|' + |7„(0|'=l, (6) 



a„(/)5:(0 + iS:(/)r„(0 = 0, a„(/)r:(0 + iS:(/)5„(/) = 0, (7) 



from which, it follows that 



|«„(Or = |iS„(/)|', \yM' = \5M\ (8) 



The "initial" conditions are 
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«„(0)=)S„(0) = 1, 7^(0) = 5„(0) = 0, 



(9) 



and )S_i(0) = 1. The transforaied Hamiltonian is given by 



^(0 = £ [A„(l)G,G_\n){n\ + B„_,(l)G_G,\n){n\ 

+ C„(l)G^ \n){n + 1| + Cl(l)G_ \n + lXn|], (10) 



where 



COr 



i«„(or-ir„(or 



+ €0 



«+|r„(of] 

+k^f^[a„(l) 71(1) + aim 



(11) 



(Or 



S„(l)\'-\P„(l)\'] + 0)[n + \P„il)\'] 



(12) 



c„(/)=^K(/)5:(/)-)s:(/)r„(/)]+«is:(/)7„(/) 

+AV^[a„(/)jS:(/) + y„(/)5:(/)], 



(13) 



for n = 0, 1, 2, • • •, and B_^(l) = -cOq II, provided that Eqs. (6) and (7) have been used. 
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Note that both A„(/) and B„(/) are real. The "initial" conditions are 



A„(0) = ^ + con, 



C„(0) = AV^. 



The diagonal part in Eq. (10) is given by 



H\l) = £ [A„(/)t7,t7»(n| + B„_,(l)(T_(T,\n){n\]. 

n=0 



Accordingly, the anti-Hermitian operator in Eq. (2) is calculated to be 



^(0 = L [MD - B„(0][C„(/)<7 Jn>(n + 1| - C:(l)(7_\n + l){n\], 

n=0 



leading to 



[77(/),//(/)] = |;[A„(/)-5„(/)]|2|c„(/)|'{(7,(7_|n)(n|-(7_(7jn + l)(n + l|) 

n=0 

-[A„(l)-B„(l)][c„(l)(r,\n){n + l\ + Cl(l)(r_\n + l){n\]}. (19) 



Therefore, Eq. (1) yields the following set of flow equations ( n = 0, 1, 2, • • •): 



dA(l) r II |2 

—^ = 2[A„(1)-B„(1)]\C„(1)\ , (20) 



dl 



= -2[A„(l)-B„(l)]\C„(l)\\ (21) 



dC(l) r n2 

"'■ = -[A„(l)-B„(l)] C„(l). (22) 



dl 



From Eq. (22), it is obvious that the off-diagonal matrix elements in fact decay as / 
grows [note that A„(/) ^ B^(l) because of the off -resonant condition in Eq. (4)]. From 
Eqs. (16) and (22) as well as the realities of A„(/) and 5„(/), C^(l) also turns out to 
be real. The exact solutions to these equations under the "initial" conditions in Eqs. 
(14)-(16) are 



A„(l) = co(n + ^] + ^L„(l), (23) 



V 



2) 2 
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5„(/) = © 



(24) 



(25) 



for n = 0, 1,2, • • where + 1) is the Rabi frequency, and L^(l) is 

given by 



M0 = 



si: 



- + 



(26) 



expi 



We note that without losing generality A„(/) is taken to be larger than BJJ). (If A 
were taken negative, then the case A„(/) < BJJ) could be realized. This point does not 



affect the subsequent discussion.) 



Furthermore, the equalities, T](l) = [dU(l) / dl]U\l) and dU(l) / dl = ri(l)U(l), 
together with Eqs. (6)-(9) and (18), lead to the fact that «„(/), jS„(/), /„(/), and 
5„(/) are all real and satisfy the relations 



«„(0 = iS„(/), 7„(/) = -5„(/), 



jS_i(/) = l. 



(27) 
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From these relations as well as Eqs. (6)-(8) and (27), the coefficients appearing in the 
unitary operator in Eq. (5) are found to be given as follows ( n = 0, 1, 2, • • •): 



a„(l)=l3„(l) = J^^, (28) 



r„(/) = -5„(/) = J— (29) 



where 



4(0 = |t4(0 + 



V "7 



(30) 



Recalling the off -resonant condition in Eq. (4), Eqs. (23)-(25) respectively converge. 



in the limit /—><», to 



( l\ ^„ 
A„(oo) = (y n + - 



(31) 



B„(oo) = co n + -]-^^E_(n), 



(32) 



C„(°o) = 0, 



(33) 
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with the energy eigenvalues, E^(n), corresponding to \E_^_j = U\oo)\e)\n) and 
Ie'.^ = U\oo)\g)\n + 1). Eq. (33) explicitly shows that the Hamiltonian is in fact exactly 
diagonalizable. 

Closing this section, we point out that although the present flow equation approach 
assumes the off -resonant case as in Eq. (4), the limiting quantities, A„(oo) and 5„(°°) 
(which are the energy eigenvalues in the states, \e)\n) and |g)|n + l), respectively) 
converge in the resonant limit to 



which are the correct results obtained by other methods [7,8]. 



IV. POSITIVE OPERATOR-VALUED MEASURES AND SUBSYSTEMS 

In this section, we discuss the effects of the unitary operation by [/(/) on the 
subsystems, i.e., the photons and atom. We shall explicitly calculate the positive 
operator-valued measures (POVMs). 

Consider |^) = |e)|n) with n^O, which is a natural "initial" state to be prepared 




(34) 




(35) 
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since in a realistic experimental situation using a cavity the photons and atom may 
initially be uncorrelated each other. The corresponding density matrix of the total 
system reads 



P^ou = \e){e\®\n){n\. (36) 



Recalling the operation of U{1) on the Hamiltonian, H(l) = U(l)HU\l), we see the 
corresponding operation on the state to be 



\^(l)) = UHl)\e)\n). (37) 



Therefore, the above density matrix transforms as 



Ptotal^Ptotal(0 = f/^(OPtota,f/(0. (38) 



The reduced density matrices of the subsystems are given by Pphoton (0 = Tr^j^^ p^^^^ (I) 
and PatomCO = Trp^^j^^P(^j^(/). The form in Eq. (36) allows us to obtain the following 
Kraus representations [10]: 



Pphoton (0= i;A,.(/)|n>(n|A;(/), (39) 



i=g,e 
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Pato»(0 = £5J0k)(e|5:(0. (40) 



m=0 



Here, {A,(/)},^^ , and {5„(/)}„.o,i,2,. fonnPOVMs 



j;a;(/)a,.(/)=v,„, (41) 



i=g,e 



Bl(l)B„(l) = I^^^^, (42) 



m=0 



with /'s being the identity operators. Using Eqs. (5) and (36), they are explicitly 
calculated as follows: 



A, (/) = {i\U\l)\e) = £ [«„(/) 5Jn){n \ + rM)SJn + l){n\], (43) 



n=0 



B„(l) = {m\U\l)\n) 



These operators do not satisfy the unital conditions: ^ A.(l) Aj (l)^ I 



photon ' 
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^ B^(l) BI (l)^ I^xom- Accordingly, entropy, which will be discussed in the next section, 

does not increase under the operations of these POVMs [12,13], in general: that is, 
whether entropy increases by the operations depends on choice of a state to be 
transformed. 

V. STATIONARY PHOTON- ATOM ENTANGLEMENT 

As mentioned in the Introduction, the flow equation approach may offer the optimal 
strategy for generating the stationary state of the Hamiltonian. This implies that it is 
possible to measure how quantum- state engineering makes a given "initial" state close 
to the stationary state. 

As mentioned in the previous section, taking into account a realistic experimental 
situation using a cavity, a natural "initial" state is a product state. Since U (I) in Eq. (5) 
is a nonlocal unitary operator, it necessarily creates entanglement between the photons 
and atom. Then, it is of interest to quantify the degree of entanglement in the engineered 
state as well as its closeness to the exact stationary state. Here, we discuss this issue by 
analyzing the entanglement entropies. 

The entanglement entropies are defined by the von Neumann entropies of the 
subsystems: 

"^photon (0 = "Trp^pjgjj P photon (0 P photon (0 > (45) 
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atom 



(46) 



where Pphoton(0 and Patom(0 are the reduced density matrices defined in the previous 
section. These are nonzero because |^(/)) in eq. (37) is an entangled state for />0. 
From the Araki-Lieb inequalities [14] in the present case. 



(47) 



it follows that the entanglement entropies of the photons and atom are identical 



"^photon (0 - 'S'atom (0 



(48) 



at each value of /, since S^^^^(l) vanishes because of purity of the total state in Eq. (38). 



Therefore, henceforth we consider only (I) . 



An explicit calculation yields 



S,,,^(l) = -s,(l)\ns,(l)-s_(mns_(l), 



(49) 



1±- 



(50) 
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where i„(/) is given in Eq. (30). Numerical analysis of S^^^^(l) with the use of the 
physical parameters taken from real experiments (for example, as in Refs. [15,16]) 
shows a monotonic behavior, exhibiting how the degree of entanglement increases as 
the "initial" product state approaches the exact stationary state. The flow of S^^^^(l) 
can be also used for quantifying the deviation of an engineered state from the stationary 
state. 

Since |^(/)) with a finite value of / is not the stationary state, its time evolution 
makes the entanglement entropies oscillatory. To see this, we have calculated the time 
evolution of 5,,,,^ (/, t) for | ^ (/, 0) = e | ^ (/)) . The result is 



Numerical analysis of S^^^^ (/, t) may show how the oscillation decays as / increases. 



■^atomC^' = t)lns^(l, t)-s_(l, t)lns_(l, t). 



(51) 




±(l-cosD.j) 




(52) 
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VI. COMMENT ON IMPLEMENTAION OF 
THE UNITARY OPERATION 

Finally, we discuss how to approximately implement the unitary operation appearing 
in the method of flow equation for the Jaynes-Cummings model. 

In a real experimental condition [15], A / A = 0(10"^). Therefore, expanding the 
coefficients, a's, fi's, y's, and d's in Eq. (5) with respect to it, we obtain the 
following approximate expression: 

U(l) = /photon ® /ato. + (l -e-'^') f « - ^-«')- (53) 

The above operator is linear in the photon field, and therefore it may not be difficult 
to be implemented in terms of the quadrature operators. On the other hand, for the atom 
operators, we note that the use of the Ramsey zone [17] yields the transformation, 
\e) ^[\e) + \g)) / ^|2 , which can be generated by the following unitary operator: 

Re,^ = ;^{^+^- + V<7, + e'V + (54) 

where the phases satisfy the condition, d-(j)=7t mod(2;r). Therefore, we find the 
atom operators to be given as follows: 
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+ V 0, ;r/ 71,0 0,71 : 



atom ' 



(55) 



(56) 



Therefore, the unitary transformation of the flow equation can approximately be 
implemented by making use of the quadrature operators of the photon field and the 
Ramsey zone. 



We have discussed stationary photon-atom entanglement in the Jaynes-Cummings 
model by applying the method of flow equation. We have explicitly constructed the 
nonlocal continuous unitary operator and the associated positive operator-valued 
measures for the photons and atom as the subsystems. We have analyzed flow of the 
entanglement entropies and furthermore have discussed an approximate implementation 
of the unitary operation. 

As emphasized in the Introduction, it is our opinion that Wegner's flow equation is 
more than just a mathematical tool for diagonalizing a given Hamiltonian. In the present 
work, we have shown how it offers a useful approach for quantum-state engineering. 
We believe that this method may be able to play a peculiar role in quantum information 
technology. 



VII. 



CONCLUSION 
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